Exact hydrodynamics of a trapped dipolar Bose-Einstein condensate 
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We derive tlie exact density profile of a iiarmonically trapped Bose-Einstein condensate (BEC) 
which has dipole-dipole interactions as well as the usual s-wave contact interaction, in the Thomas- 
Fermi limit. Remarkably, despite the non-local anisotropic nature of the dipolar interaction, the 
density turns out to be an inverted parabola, just as in the pure s-wave case, but with a modified 
aspect ratio. The "scaling" solution approach of Kagan, Surkov, and Shlyapnikov [Phys. Rev. 
A 54, 1753 (1996)] and Castin and Bum [Phys. Rev. Lett. 77, 5315 (1996)] for a BEC in a 
time-dependent trap can therefore be applied to a dipolar BEC, and we use it to obtain the exact 
monopole and quadrupole shape oscillation frequencies. 
PACS numbers: 03.75.Kk, 34.20. Cf, 32.80. Qk, 75.80.+q 



The current revolution in the cooling of atomic gases to 
ultra-cold temperatures has led to the experimental re- 
alization of a number of highly idealized and fundamen- 
tal regimes in the quantum physics of many-particle sys- 
tems ^ 12 ■ For example, in a Bose-Einstein condensed 
gas containing several million atoms the zero-point kinetic 
energy can be negligible in comparison to the interpar- 
ticle interaction energy. In this so-called Thomas-Fermi 
regime the collective dynamics of the Bose-Einstein con- 
densate (BEC) are described by the coUisionless hydrody- 
namic theory of Bose superfluids at zero temperature |3] . 
Furthermore, the interparticle interactions can be directly 
manipulated using external electromagnetic fields, giving 
a degree of freedom not normally available in the systems 
traditionally studied in statistical physics. The dominant 
interactions in these gases are usually of the van der Waals 
type, decaying as r~^, short-ranged in comparison to the 
interatomic separation. This permits an effective descrip- 
tion of the interactions solely in terms of the s-wave scat- 
tering length, Og. Within the mean- field scheme of the 
Gross-Pitaevskii equation ^ the interactions are modelled 
by the contact pseudo-potential, g5(r) = (47ras?i^/m)(5(r), 
where m is the atomic mass, and g incorporates the quan- 
tum aspects of low-energy scattering. By using a mag- 
netic field to induce a Feshbach resonance during atomic 
collisions, one can adjust the value of Og between positive 
(repulsive) and negative (attractive) values [1]. 

Here we consider a harmonically trapped BEC which, as 
well as van der Waals interactions, also has dipole-dipole 
interactions. Such a system was first considered by Yi 
and You [H] in the context of dipoles induced by an ex- 
ternal electric field and by Coral, Rzazewski, and Pfau 
[7] in the context of permanent magnetic dipoles aligned 
by an external magnetic field. In comparison to van der 
Waals interactions, dipole-dipole interactions have a much 



longer range and are anisotropic. This has been pre- 
dicted to lead to BECs with unusual stability properties 
IHlini, exotic ground states such as supersolid j^E] and 
checkerboard phases and modified excitation spec- 
tra ^1 even to the extent of a roton minimum |2| ■ 
Dipole-dipole interactions are also inherently controllable, 
either via the magnitude of the external electric field in the 
case of electrically induced dipoles, or in the case of per- 
manent dipoles by modulating the external aligning field 
in time ^Hl) which allows control over both the magnitude 
and sign of the dipolar interactions. 

The long-range part of the interaction between two 
dipoles separated by r, and aligned by an external field 
along a unit vector e, is given by 
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where the coupling Cdd/(47r) depends on the specific real- 
ization. Electric dipoles induced by an external static elec- 
tric field E = Ee have Cdd = E'^o? je^ jS], where a is the 
static polarizability. For magnetic atoms with magnetic 
dipolc moment dm aligned by an external magnetic field 
B = Be, then Cdd = t^-od^i- A measure of the strength 
of the dipole-dipole interaction relative to the short-range 
van der Waals s-wave scattering energy is given by the 
dimensionless quantity 



Edd 



Cdd 



(2) 



For edd > 1 instabilities are expected [TIE] in the limit of 
negligible kinetic energy. Note that the value of g can be 
altered in the presence of a strong external electric field 
|12| . Alkali atoms can have a magnetic dipole moment 
of dm = (Bohr magneton), and for ^"^Rb, which has 
as « 103ao (Bohr radius), one finds Edd ~ 0.007. Na has 



£dd ~ 0.004, and accordingly magnetic dipolar effects in 
BECs of these atoms are expected to be small, at least in 
a stationary condition. The effects of the dipolar interac- 
tions can be made visible by rotating the magnetic field 
in resonance with a collective excitation frequency of the 
system 15 . Also, using a Feshbach resonance to reduce 
Og could substantially increase edd |lfc>j . However, a good 
candidate atom for a dipolar BEC is chromium which has 
dm = 6/iB- The s-wave scattering length of ^^Cr, the 
more common isotope, is 170 ± 40 ap. Consequently for 
that isotope Cdd ~ 0.089. The less common isotope ^°Cr 
has fls = 40 ± 15 ao, and £dd ~ 0.36 is much higher. A 
chromium BEC does not yet exist, but a program of ex- 
periments is underway with that goal in mind |17| . Using 
a crossed optical trap, like in the recent cesium condensa- 
tion experiment 18 , atom numbers in a ^^Cr BEC on the 
order of iV = 10^ might be within reach, and the radial 
and longitudinal frequencies of the cylindrically symmet- 
ric (optical) harmonic confinement would be in the region 
oiujx — ujy = 27rl70 s^^ and ujz — 27r240 s^^, respectively 
|17| . The trap anisotropy 7 = w^/wj, = 1.41 is optimal for 
enhancing the condensate shape deformations ^Hl- For a 
harmonically trapped BEC with just s-wave scattering the 
Thomas-Fermi limit is reached for large values of the di- 
mensionless parameter iVos/aho |2] where Oho — y^h/muj 
is the harmonic oscillator length of the trap, and N is the 
total number of atoms. For ^^Cr one has Nas/aho ~ 90, 
illustrating that the Thomas- Fermi regime is an appropri- 
ate starting point for analyzing these experiments. 

The application of coUisionless hydrodynamics to an or- 
dinary BEC is based upon the existence of a macroscopic 
condensate order parameter ■(/'(r, t) = \/n{r, t) exp i0(r, t), 
normalized to consisting of density n(r,t), and phase 
(/)(r, t) . The phase determines the superfluid velocity via 
potential flow, v(r,t) — {h / m)W 4){r , t) . The nonlinear 
evolution of the condensate in the Thomas-Fermi regime 
is then dictated by the continuity and Euler equations 
given by, respectively, 

-V.(nv), (3) 



dt 
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Vc-x_t is the external potential, which is provided by the 
harmonic trap, ycxt(r) — m{ujlx^ + ujyy"^ + cj^z^)/2. Eqs. 
(O and ^ describe the potential flow of a fluid whose 
pressure and density are related by the equation of state 
P = (.g/2)n^, and are equivalent to the time dependent 
Gross-Pitaevskii equation |4| when the zero-point fluctu- 
ation term {—h^V'^y/n/2rn^/ri) is included in the brack- 
ets on the right hand side of Q after the Bernoulli term 
TOW^/2. The equilibrium solution of Q and Q), which has 
V = 0, is the well known Thomas- Fermi inverted parabola 
density profile 

(5) 



n{v) = - V;xt(r)]/.9 for n{v) > 0, 



and n(r) =0 elsewhere. Thus the density profile is com- 
pletely determined by the trapping potential, having the 



same aspect ratio as the trap. 

An analysis of the dynamical (v 7^ 0) solutions of © 
and (0J) has been given by Stringari |3j. However, mo- 
tivated by the experimental possibility of changing the 



trapping frequencies in time, so that ujj 



j(t) with 



J — x,y,z, Kagan, Surkov, and Shlyapnikov 27 and 
Castin and Dum studied a special class of "scaling" 
solution, valid for time-dependent harmonic traps, corre- 
sponding to [2] 
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v(r,t) = [aa^{t)x^ + ay{t)y'^ + az{t)z^ 



, (6) 
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valid where n(r, i) > and n{r,t) = elsewhere. The 
central density no{t) is constrained by normalization to 
be no{t) = 15N/[8TTRx{t)Ry{t)R,{t)]. The time-evolution 
of the radii Rj reduces to the solution of three ordi- 
nary differential equations, and the components of the 
velocity field, aj, depend on the corresponding radii as 
aj = Rj/Rj. Although these scaling solutions represent 
only a limited class of the possible solutions to the hydro- 
dynamic equations, they model two crucial experimental 
situations: i) the expansion of a BEC when the trap is re- 
laxed or turned off — this is usually necessary to image the 
BEC which would otherwise be too small in situ (time of 
flight measurements then record the original momentum 
distribution); ii) the response of the BEC to a modula- 
tion of the trapping frequencies, which is the easiest way 
to generate shape oscillations of the BEC. 

The existence of the class © of analytic scaling solu- 
tions relies upon a harmonic trap and the local character of 
s-wave interactions. The remainder of this paper demon- 
strates an analogous exact result for a dipolar BEC. This 
result is not obvious because of the non-local nature of 
dipolar interactions. Nevertheless, it turns out that the 
dipolar contribution to the mean-field potential generated 
by a parabolic density profile is itself also parabolic so that 
(jnj and still represent a self-consistent solution to the 
hydrodynamic equations of motion. 

Exact static solution for a Thomas-Fermi dipolar BEC 
— Dipolar interactions can be included in a mean-field 
treatment of the BEC via a Hartree potential 



$dd(r) 
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as shown by Yi and You jj^l based on a calculation of the 
two-body T-matrix. The hydrodynamic equations of mo- 
tion for dipolar BEC are then the same as above but with 
the Euler equation Q supplemented by adding <i>dd(r) 
into the bracket on the r.h.s. We begin our analysis with 
a suggestive re-casting of the dipole-dipole term 
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We can then write 



$dd(r) = -Cdd e,ejV,Vj(^(r) + -n{r) (10) 



with 



d^r' n(r') 



(11) 



The problem thereby reduces to an analogy with electro- 
statics, and one need only calculate the 'potential' 0(r) 
arising from the 'static charge' distribution n{r). We 
are now in a position to see the basic form n(r) must 
take, since <^(r) given by obeys Poisson's equation 
= — n(r). Thus, if the density is parabolic, of the 
form |(BJ, then Poisson's equation is satisfied by a poten- 
tial of the form (p = gq + aix^ -\- a2y^ + a^z^ + 04x^2/^ -f- 
a^x^z^ + ae.y'^z'^+ajx'^ + asy'^ + agz'^, and by ((Tn|l $dd(r) is 
also parabolic. Therefore, an inverted parabola remains a 
self-consistent solution to the dipolar hydrodynamic equa- 
tions. The evaluation of the integral Hll() is too long to be 
reproduced here |22|, but can be achieved either by inte- 
grating over successive thin ellipsoidal shells, or in the case 
of cylindrical symmetry by transforming into spheroidal 
coordinates and using the known Green's function for Pois- 
son's equation in that basis. 

For simplicity (though not necessity), we shall restrict 
attention to cylindrical symmetry for both trap and BEG 
about the z axis, along which the external polarizing field 
is understood to point. Thus we have Rx — Ry for the 
BEG and ^oxt = [■m/2)[ujlp^ + ujIz\ where ^ + 
y^, for the trap. Then, inside the condensate region, the 
dipolar potential (|10() generated by an inverted parabola 
density distribution ^ of dipoles, is [221 
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with K = Rx/Rz being the condensate aspect ratio. We 
recover the result (9) of in the particular case of a 
spherical dipolar condensate {Rx = Rz)- In order to de- 
termine the equilibrium values of Rx and Rz, the result 
((T^ for $dd should be substituted into the Euler equation 
Q with V = 0, giving a transcendental equation for the 
aspect ratio k 



3k Edd 



^ + 1 



/(«) 



- 1 + (£dd - 1) («' - 7') = 

(14) 

where 7 = Uz/ujx is the ratio of the harmonic trapping fre- 
quencies. Fig. ^ shows some examples of the dependence 
of K upon edd for oblate, spherical and prolate traps. The 
effect of the dipole-dipole forces polarized along the z-axis 
is to make the condensate more cigar-shaped along z. For 
an oblate trap (7 > 1) the BEG becomes exactly spher- 
ical when edd = (5/2) (7^ — l)/(7^ + 2). An analysis of 
the energy functional associated with the scaling solution 
(0 shows that when Sdd > 1 the solution loses its global 




Figure 1: Equilibrium aspect ratio k of the BEG, as a 
function of the dipole-dipole to s-wave coupling ratio, £dd- 
Each line is for a trap of a different aspect ratio 7 (k = 7 
when £dd = 0). When < 7 [k] < 1 the trap [condensate] 
is prolate, 7 [k] > 1 the trap [condensate] is oblate. The 
dashed lines indicate an unstable branch. However, local 
density perturbations most likely render the condensate 
unstable whenever ^dd > 1- 



stability to scaling perturbations and is only metastable 
8 , 12 , 22 (only a local minimum of the energy), the global 
minimum being a collapsed pencil-like prolate state. In- 
creasing Edd yet further, the scaling solution ceases to exist 
entirely in the region marked as unstable in Fig. ^ How- 
ever, we strongly caution the reader that any state of a 
dipolar BEG with Edd > 1 is probably unattainable since 
other classes of perturbations (that do not respect the scal- 
ing form ©), such as local phonons, render the dipolar 
BEG unstable to collapse whenever £dd > 1 [71 ITS). 

A moment's thought shows that the aspect ratio is in- 
dependent of the specifics of the density profile and indeed 
the same transcendental Eq. H14I) has been previously ob- 
tained using a Gaussian ansatz |12| . However, what is 
specific to the exact solution are the absolute radii which, 
once ()14|) has been solved for k, are given by 
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and Rz — Rx/k- If the trap is turned off the s-wave 
and dipole-dipole interaction energies are converted into 
kinetic energy, the so-called release energy i^rci, which can 
be directly measured in an experiment. One finds = 
15gN^[l - eMfii^)]/i28TrRlRz). 

Exact scaling hydrodynamics of a dipolar condensate — 
The parabolic form of $dd(r) means that the scaling so- 
lutions |21)[I2H for dynamics also apply to a dipolar BEG. 
Substituting the scaling solutions JB)) and Q into the con- 
tinuity and Euler equations yields the following classical 



equations of motion for the radii 
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Solving these equations under particular choices of the de- 
pendence of LOj{t) upon time gives the evolution of the 
density and phase of the condensate when undergoing e.g. 
free expansion (trap switched off), or shape oscillations 
(time-modulation of the trap). A time-dependent dipo- 
lar coupling is also possible as explicitly indicated. As 
mentioned above, our results show that previous calcula- 
tions based upon the Gaussian variational ansatz of the 
aspect ratio K{t) during ballistic expansion |19j . or during 
the parametric excitation of the quadrupole shape oscilla- 
tion mode using dipolar interactions as a time-dependent 
coupling |15| . are exact in the Thomas- Fermi limit. 

Linearizing UlSHlTfl around the equilibrium solution 
(|15|l gives the following frequencies for small amplitude 
scaling oscillations 
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ri- corresponds to quadrupole, and fl+ to monopole 
(breathing) shape oscillations, respectively EIj see 
Fig. 121 Note that the frequencies do not depend on g 
directly, only through Edd- Fig. |31 illustrates that certain 
trap anisotropics maximise the effects of the dipolar inter- 
actions upon the collective excitation frequencies. 

Conclusion — Dipolar condensates promise to open a 
new chapter in the study of controllable superfluids. At 
first sight, the non-local and anisotropic form of dipolar 
interactions also make the elucidation of the static and 
excited states of a dipolar BEG non-trivial, i.e. direct nu- 
merical calculation is for instance computationally heavy 
in comparison to standard BECs with contact interactions. 
We have shown, however, that there exist simple exact so- 
lutions which should prove useful for quantitative analysis 
of upcoming experiments. 
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Figure 2: The monopole (dashed) and quadrupole (dash- 
dotted) oscillation frequencies in units of the radial trap- 
ping frequency, as functions of e^d- The solid lines are for 
a BEG with only s-wave contact interactions, il^.. The 
trap aspect ratio is set at 7 = 240/170. The frequencies 
become complex at precisely the value of e^d a-t which the 
equilibrium solution becomes unstable. However, other in- 
stabilities, due to phonons, are possible whenever e^id > 1- 
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Figure 3: The monopolar (dashed) and quadrupolar 
(solid) frequencies as fractions of their s-wave-only fre- 
quencies, nl, of 52Gr (edd = 0.089) and ^"Cr (edd = 0.36), 
plotted as functions of the trap anisotropy 7. 
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